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ABSTRACT 


This  report  constitutes  the  final  report  of  research  with  ONR  under 
contract  no.  N00014-79-C-0555,  and  entitled  "Conversion  Efficiencies  of 
Charged  Particle  Beam  Driven  Free  Electron  Lasers". 

The  contract  covered  the  period  from  June  11,  1979  to  November  30, 

1984.  Technical  progress  herein  reported  covers  this  entire  period;  however, 
in  depth  technical  details  specialize  in  the  final  year  as  other  annual 
reports  have  thoroughly  reviewed  the  work  of  previous  years.  Additionally, 
the  totality  of  publications  and  talks  resulting  directly  from  this  contract 
over  the  entire  contract  period  are  listed. 
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I.  INTRODUCTION 

The  general  purpose  of  the  research  under  this  contract  was  to  develop 
sophisticated  models  of  electron  beam  driven  free  electron  lasers  in  order 
to  more  accurately  predict  radiation  gain  rates  and  fundamental  efficiencies. 
Over  the  period  covered  by  the  contract  a  number  of  magnetic  field  and  beam 
configurations  were  studied  but  primary  emphasis  has  been  on  helical,  longi¬ 
tudinal,  and  planar  wigglers  with  strong  guide  fields. 

Within  this  area  low  energy,  high  current  devices  were  studied  the  most, 
aiming  at  an  understanding  of  high  gain,  long  pulse  devices.  Such  systems 
usually  formulate  the  analytic  amplification  problem  by  choosing  beam  equilib¬ 
rium  and  linearly  perturbing  about  it  to  calculate  (linear)  gain  rates,  and 
then  the  nonlinear  saturation  mechanisms  are  studied  separately  in  order  to 
evaluate  the  leading  limiting  mechanism.  Again  such  analyses  dominated  our 
research  program  and  these  results  are  discussed  in  the  following  sections, 
previous  progress  reports,  and  our  total  publications.  Additionally,  part 
of  our  studies  analyzed  the  all-important  question  of  the  existence  of  an 
equilibrium  about  which  to  perturb  and  significantly  it  was  found  that  the 
properly  posed  problem  of  a  beam  with  finite  emittance  being  emitted  from 
a  -diode  does  not  always  give  a  satisfactory  equilibrium  type  solution. 

This  example  has  pointed  out  the  need  for  more  careful  analyses  of  the  so- 
called  equilibrium  problem.  As  FEL  gain  rates  tend  to  be  sensitively  depen¬ 
dent  on  the  shape  of  the  velocity  space  distribution  functions  and  nonlinear 
saturation  depends  on  the  three  dimensional  nature  of  the  amplified  wave, 
wiggler,  and  beam  distribution  function  a  properly  chosen  set  of  beam  propa¬ 
gation  features  must  be  carefully  determined.  Most  probably  this  will 
require  a  very  strong  coupling  between  relevant  theory  and  experiment. 
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Finally  it  should  be  noted  that  the  research  conducted  under  this 
contract  laid  the  foundations  for  research  on  other  FEL  configurations, 
namely  those  of  the  high  energy,  zero  guide  field  problems  of  LANL  and  LLNL 
and  magnetic  crystal  FEL  Y  radiation  generation. 

In  the  remainder  of  this  report  details  of  the  supported  research 
over  the  last  contract  year  are  presented  in  Section  II, and  Section  III 
lists  contract  supported  publications  and  oral  presentations. 
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II.  RECENT  PROGRESS 

In  this  section  we  review  the  progress  made  under  this  contract  over  the 
last  year. 

A.  Profile  of  a  Relativistic  Electron  Beam  Propagated  through  a  Linear 

Wiggler  and  an  Axial  Guide  Field 

Most  theoretical  analyses  of  the  free-electron  laser  instability  employ 
highly  simplified  models  in  which  the  electron  density  is  taken  to  be  uniform 
in  space.  Noting  that  a  function  of  the  constants  of  motion  of  a  single 
electron  is  a  steady-state  solution  of  the  Vlasov  equation,  Davidson  and  Uhm1 
were  the  first  to  construct  self-consistent  equilibria  for  electron  beams 
propagated  in  an  axial  guide  field  and  a  helical  wiggler  field,  and  the 
resulting  spatial  modulation  of  the  electron  density  was  determined  for 
different  choices  of  the  functional  dependence  of  the  distribution  function 
on  the  constants  of  motion. 

In  this  work,  we  have  examined  the  experimental  situation  where  the 
directional  flux  of  electrons  (and  beam  emittance)  at  the  anode  is  specified 
and  the  beam  is  allowed  to  propagate  through  a  drift  tube  immersed  in  a  strong 

axial  guide  field  and  a  linear  (or  planar)  wiggler.  Using  the  nonlinear  orbits 

2 

obtained  earlier  at  SAIC/PRI,  the  distribution  function  of  the  electrons  at  any 
point  along  the  drift  tube  is  determined  from  the  equation  of  continuity  in 
phase  space  and  the  electron  density  is  then  given  by  an  integration  over 

2 

momenta.  Unlike  the  calculation  in  Ref.  1,  the  orbits  determined  earlier 
permit  us  to  treat  the  case  where  there  is  a  resonance  between  the  gyromotion 
of  the  electrons  and  the  effective  axial  periodicity  induced  by  the  wiggler. 
Additionally,  the  radial  dependence  of  the  wiggler  field  is  allowed  for. 

In  Ref.  2  it  is  argued  that  in  the  presence  of  an  axial  guide  field  it  is 
natural  to  introduce  the  gyroradius  p,  gyroangle  0,  and  guiding  center  coordinates 
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X  and  Y.  It  is  then  shown  that  for  sufficiently  small  wiggler  fields  (measured 
by  e  =  6H/Hq,  <5H  being  the  wiggler  and  Hq  the  guide  field  intensity)  the  motion 
of  the  electrons  can  be  locked  onto  a  single  resonance,  with  the  particles 
strongly  influenced  by  this  resonance  but  mildly  so  by  adjacent  harmonics. 

This  resonance  condition  has  the  form  SLB  +  kz  -  0  for  some  integer  l,  where 
k  =  2tt/A,  A  is  the  wiggler  wavelength,  and  z  is  the  axial  velocity.  It  is 
further  shown  that  in  this  case  the  motion  is  nonstochastic, and  explicit 
integration  of  the  equations  of  motion  is  then  possible.  For  FEL  applica¬ 
tions  it  is  imperative  to  limit  the  strength  of  the  wiggler  so  as  to  avoid 
stochastic  motion. 

3 

In  the  experiment  carried  out  by  Roberson  et  al.  at  the  Naval  Research 
Laboratory,  which  provided  the  original  motivation  for  our  work,  a  current 
of  electrons  (=*500  A)  is  accelerated  through  a  potential  difference  of 
350  kV  (y  sx  1.68)  before  impinging  on  an  anode  of  radius  R^  =*  1/2  cm  and 
entering  a  drift  tube  of  length  60  cm.  In  this  experiment  electric  and 

magnetic  self-field  effects  on  particle  orbits,  which  are  measured  by  the 
dimensionless  parameter  (2yc0pA^) [1  -  {z/c)\  -  1/2,  are  nonnegligible. 

(Here  co  is  the  plasma  frequency  and  the  guide  field  cyclotron  frequency.) 

Further,  in  the  experiment  the  magnetic  field  was  allowed  to  build  to  its 
full  value  in  a  transition  region  preceding  the  drift  tube  entrance.  The 
work  described  in  Ref.  2  and  herein  is  idealized  in  the  sense  that  neither 
self-field  effects  nor  the  effect  of  a  transition  region  is  allowed  for. 

The  electron  density  at  r  is  given  by  n(r)  =  J  d  jd  f(r,p),  where  f, 
the  electron  distribution  function  at  (£,£),  satisfies  the  continuity  equa¬ 
tion  in  6-dimensional  phase  space:  df/dt  =  0.  Knowing  the  beam  emittance 
at  the  anode  and  the  particle  orbits,  the  indicated  integration  over  momenta 
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can  then  be  performed  to  determine  the  density  as  a  function  of  distance  along 
the  drift-tube.  A  reasonable  form  for  the  beam  distribution  at  the  anode, 
located  at  z  =  0,  is 

f  =  #0(Ro  -  R)6[(y2  -  l)m2c2  -  i)  exp(-<j>)  , 


wherein  is  a  normalization  constant,  0(*)  is  the  Heaviside  unit  step 

function,  RQ  is  the  anode  radius,  R  =  (x2+y2)^  is  the  radial  distance  of 

the  point  (x,y)  from  the  axis  of  the  drift-tube,  and 

2 

_  t*  v  n  -nr 


<t>  = 


p- +  pi  t 


-2 

Pi 


Pz  -  Pb 


Pm  and  pA  being  the  spread  in  the  parallel  and  perpendicul ar  momenta  of 

the  electrons  emerging  from  the  anode,  with  the  parallel  momenta  being 

centered  about  the  resonant  value  pfa  =  j>m|Qo|/k  (Ref.  2).  The  Heaviside 

function  in  the  expression  for  f  indicates  that  a  beam  of  radius  R  and 

o 

of  uniform  density  emerges  from  the  anode.  The  6-function  indicates  that 
the  beam  is  monoenergetic,  with  a  fixed  y.  The  experimentally  determined 
spread  in  the  perpendicular  speed  is  roughly  given  by  -(pA/mc) a  0.08-0.24;  by 

p 

a  simple  kinematical  argument  one  finds  that  p„  =  %(pA/p  ),  whence  the 
spread  in  parallel  speed  is  determined. 

In  Fig.  1  the  beam  density  (in  arbitrary  units)  is  plotted  as  a  function 
of  the  axial  distance  z  from  the  anode  along  the  drift  tube  for  e  =  1/1024. 

For  the  upper  curve  the  observation  point  at  which  the  density  is  measured  has 
coordinates  x  =  0,  y  =  0.5  while  the  lower  curve  has  x  =  0,  y  =1.0x10  .  All 

linear  dimensions  are  given  in  centimeters  and  plotted  results  are  accurate 
to  <  1%.  The  substantial  density  modulation  evident  in  the  upper  curve  can 
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be  understood  by  examining  the  nature  of  the  particle  orbits  using  Fig.  2. 

In  this  diagram  circle  A  represents  the  anode  (radius  1.2)  and  0  is  the 
observation  point  (coordinates  x  =  0,  y  =  0.5).  Consider  gyro-orbits  which 
enclose  circle  R,  pass  through  observation  point  0  and  have  their  guiding 
center  on  the  y-axis  (circle  R  for  example).  Particles  in  such  orbits 

X/ 

necessarily  move  beyond  the  anode  radius  during  part  of  their  trajectory. 
Since  the  beam  has  little  spread  in  parallel  momentum  by  far  the  greatest 
number  of  particles  have  gyroradii  of  0.453,  that  of  particles  exactly  at 
resonance  (p  =  p  =  £m[ft  |/k).  Circles  R  and  R  are  gyro-orbits  that 
bound  this  most  probable  value,  and  are  populated  with  1  x  10-^  of  the  number 

of  particles  at  resonance.  As  we  move  forward  in  space  from  the  anode,  the 

particles  in  the  relevant  orbits  (denoted  by  the  shaded  regions)  move  outside 
A  and  the  density  at  0  falls.  For  £  =  2,  the  resonance  condition  £§  +  kz  =  0 
implies  that  in  one  wiggler  wavelength  (3  cm)  the  particles  complete  about 
half  of  their  gyro-orbit  and  the  density  should  be  minimum.  Over  a  distance 
of  two  wiggler  wavelengths  the  particles  return  to  their  original  locations 
(in  the  limit  where  all  parallel  momenta  are  identical)  and  the  density  is 
again  at  its  peak. 

Because  a  small  spread  in  parallel  momenta  actually  exists,  particles 
completing  their  orbits  return  to  their  original  positions  at  slightly  dif¬ 
ferent  times,  producing  a  smaller  modulation  of  density  evident  in  the 
shorter  wavelength  oscillations  in  the  upper  curve  of  Fig.  1  and  the  latter 

part  of  the  lower  curve.  (The  initial  osci 1 1 ations  of  density  in  the  lower 

curve  are  numerical  noise.) 

Of  course,  there  are  other  regions  besides  the  shaded  area  in  Fig.  2 
which  contribute  to  the  plotted  density  modulations.  We  have  simplified 
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our  discussion  and  diagram  to  keep  the  argument  clear.  By  recalling,  however, 
that  the  only  significant  gyroradii  are  approximately  0.45  and  that  the  obser¬ 
vation  point  is  0,  it  is  easy  to  see  that  the  only  particles  contributing  to 
the  density  modulation  lie  in  an  appropriately  expanded  shaded  region  symmet¬ 
rically  placed  about  the  y-axis.  For  instance,  since  the  radius  of  circle  L 
in  Fig.  2  is  0.85,  practically  no  particle  with  a  guiding  center  lying  on  the 
y-axis  to  the  left  of  0  whose  orbit  passes  through  0  moves  out  of  A.  Such 
a  particle,  therefore,  does  not  contribute  to  the  observed  density  depletion. 

The  wiggler  acts  to  alter  the  speed  of  the  particles  as  they  move  along 
the  drift  tube.  By  increasing  the  amplitude  of  the  wiggler  field  an  enhance¬ 
ment  of  the  shorter-wavelength  modulation  seen  in  Fig.  1  is  expected.  This 
is  clearly  evident  in  Figs.  3  and  4  in  which  e  =  1/128  and  e  =  1/64  respectively. 
It  is  this  short-wavelength  modulation  that  is  apparently  discussed  in  Ref.  1, 
for  in  that  work  it  is  assumed  there  is  no  spread  in  the  parallel  momenta  of 
the  electrons. 

We  next  define  the  parallel  temperature  by 


nTz  e 


d  £  f (r,£)ym 


p  2 

—  -  V  (r) 
ym  z | 


where  the  mean  parallel  velocity  is  given  by 


nVz  = 


d  £ 


The  quantity  Tz  is  a  measure  of  the  spread  in  parallel  velocities  in  the  beam 

frame.  In  order  that  the  electromagnetic  radiation  emitted  by  the  electrons 

be  as  coherent  as  possible,  Tz  must  be  small.  The  normalized  parallel  temperature 
0  9 

T  =  (ck/S,|^Q|)  (yTz/mc  )  as  a  function  of  z  for  e  =  1/128  and  y  =  0.5  is  shown 
in  Fig.  5.  The  initial  rise  in  T  can  be  understood  as  the  consequence  of  an 
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effective  spread  in  parallel  momenta  as  particles  with  pz  -  p^  orbit  out  of 
the  anode  radius  and  vacate  observation  point  0.  The  subsequent  increase  in 
T  is  due  to  asynchronism  of  particle  motion.  Note  that,  although  T  increases 
by  almost  two  orders  of  magnitude  it  is  still  not  very  large,  about  10  keV 
at  z  -  24  cm. 

The  effects  of  increasing  the  beam  energy  on  the  density  modulation 
are  shown  in  Fig.  6.  Here  e  =  1/128,  y  =  0.5  and  the  %  = 3  harmonic  is 
resonant  (y  =  2.62).  Note  the  larger  density  depletion  than  in  the  1  =  2 
case.  This  is  readily  understood  since  the  most  probable  gyroradius  is 
now  0.65  and  a  greater  fraction  of  electrons  have  trajectories  which  take 
them  beyond  the  anode  radius. 

In  conclusion,  this  work  indicates  that  significant  density  modulation 
can  be  expected  from  the  axial  guide  field,  the  wiggler,  and  the  spread  in 
parallel  momenta.  This  may  be  expected  to  greatly  limit  laser  gain,  par¬ 
ticularly  in  cases  where  a  typical  gyroradius  is  not  appreciably  smaller 
than  the  anode  dimension. 

B.  Influence  of  Finite  Radial  Geometry  on  Coherent  Radiation  Generation 

by  a  Relativistic  Electron  Beam  in  a  Longitudinal  Magnetic  Higgler 

The  influence  of  finite  radial  geometry  on  the  longitudinal  wiggler  free 
electron  laser  instability  has  been  investigated^  for  TE  mode  perturbations 
about  a  uniform  density  electron  beam  with  radius  (see  Fig.  1).  The  equi¬ 
librium  and  stability  analysis  is  carried  out  for  a  think  tenuous  electron 
beam  propagating  down  the  axis  of  a  multiple-mirror  (undulator)  magnetic  field 
Bq(x)  =  Bq[1  +  (<SB/B0)sinkQz]ez,  where  =  2ir/k0 


=  const,  is  the  wiggler 
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2^ 

wavelength.  It  is  assumed  that  kgRj^  «  1,  and  that  perturbations  are  about  the 

A  n  A 

self-consistent  Vlasov  equilibrium  f^(x,p)  =  ( n^/u ) 6 ( pf  - 

2  2  2 

x  5(PZ  -  Y^V^) ,  where  Pi  =  Pr  +  p0 ,  Pf  is  the  canonical  angular  momentum,  and 
n^,  Yf-j »  ^5’  ^i^  and  are  positive  constants.  For  6B/Bg  «  1  and  slow  beam 
rotation  (co^  «  =  eBg/y^mc),  the  equilibrium  density  is  uniform  (rig)  out  to 

A  A  1/ 

the  beam  radius  Rg  =  (ZT^^/y bnftobaiCb) 2'  detailed  free  electron  laser  stability 
properties  are  investigated  for  the  case  where  the  amplifying  radiation  field  has 

^  A  /\ 

TE-mode  polarization  with  perturbed  field  components  (6EQ,  6B  ,  6B  ).  The 

0  r  z 

3 

matrix  dispersion  equation  is  analyzed  in  the  diagonal  approximation,  and  it  is 

/\ 

shown  that  the  positioning  of  the  beam  radius  (Rg)  relative  to  the  conducting 
wall  radius  (R  )  can  have  a  large  influence  on  the  growth  rate  and  detailed 
stability  properties.  Analytic  and  numerical  studies  show  that  the  growth  rate 

/s  3 

increases  as  R./R  is  increased, 
b  c 

C .  Compton  and  Raman  Free  Electron  Laser  Stab ility  Properties  for  a 
Cold  Electron  Beam  Propagating  through  a  Helical  Magnetic  wiggler 
Under  the  auspices  of  the  present  contract  we  have  given  an  extensive  charac¬ 
terization^  of  the  range  of  validity  of  the  Compton  and  Raman  approximations  to 
the  exact  free  electron  laser  dispersion  relation  for  a  cold,  relativistic 
electron  beam  propagating  through  a  constant-amplitude  helical  wiggler  magnetic 
field  Bw  =  -BcoskgZex  -  Bsink^ze^.  Here,  Xg  =  Z-rr/k^  is  the  wiggler  wavelength, 

B  the  wiggler  amplitude  (assumed  constant),  and  the  electron  beam  is  treated 
as  infinite  in  transverse  extent.  For  example,  a  detailed  numerical  analysis 
shows  that  the  Compton  approximation  (<$<J>  -  0)  gives  a  valid  estimate,  to  within 
ten  percent,  of  the  maximum  growth  rate  of  the  upshifted  emission  peak  for 
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system  parameters  satisfying 


ojpck0 


< 


Yb  U  +  eb) 
25 


U) 


2  a  2 

Here,  co^  =  4TTn^e  /ym  is  the  relativistic  plasma  frequency-squared,  ooc  =  eB/ymc 

2  -J' 

is  the  relativistic  cyclotron  frequency,  y^  =  (1  -  Bb)  2  is  the  relativistic 
mass  factor,  B^c  =  p^/ym  is  the  axial  beam  velocity,  y  is  defined  by 

o  •)  0  0  O  A  0  1 

y  =  (1  +  pQm  c  +  e  B  /m  c  kg) 2,  and  Gq(pz)  =  <S(pz  -  Pg)  is  the  equilibrium 
axial  momentum  distribution  of  the  beam  electrons. 


D .  Influence  of  Beam  Quality  on  Free  Electron  Laser  Growth  Rate  for  an  Intense 

Electron  Beam  Propagating  through  a  Helical  Magnetic  Wiggler 
We  have  made  use  of  the  fully  kinetic  dispersion  relation  derived  by 
Davidson  and  Uhm  to  investigate  the  influence  of  beam  quality  on  FEL  growth 
rate  for  a  relativistic  electron  beam  propagating  through  a  helical  wiggler 

/N  /\  ZN 

field  B  =  -BcosKnze„  -  Bsinknze  .  Assuming  a  sufficiently  tenuous  beam  that 
the  Compton-regime  approximation  is  valid  (ficf>  =  Oj,  and  assuming  weak  resonant 
instability  with  small  temporal  growth  rate  y  satisfying  |y  /&>.  |  «  1  and 

K  K  K 

lyk/kAvzl  «  1,  it  is  found  that  the  growth  rate  can  be  approximated  by 


Yk 


IT 

4 


CO 


c  2-322 
2~J  udY  m  c 


c  k 


0 


(2) 


for  small  fractional  momentum  spread  in  the  beam  electrons  (Apz/pg  «  1).  Here, 
Gg(pz)  is  the  equilibrium  distribution,  coc  =  eB/ymc  is  the  cyclotron  frequency, 
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ojp  -  ^irn^e^/ym  is  the  plasma  frequency-squared,  and  ymc2  =  (m2c4  +  c2pg  +  e2B2/kg 'yL 
is  the  average  energy. 


For  waves  excited  with  positive  phase  velocity  k  >  0,  it  follows  from 


Eq.  (2)  that  Yk  <  0  accordingly  as  [y8Gg/3p  ] 


0  HzJ  v 2  =oj k / K 


<  0.  That  is,  waves  with 


phase  velocity  in  the  region  of  positive  momentum  slope  of  Gy(Pz)  are  amplified, 
corresponding  to  instability  with  yk  >0.  As  a  specific  example,  we  have 
investigated  detailed  growth  rate  properties  for  the  case  where  the  beam  distri¬ 


bution  is  Gaussian  in  p  with 


Go<y  -  <!>\i 


exp 


2, 


Ap. 


(3) 


Typical  results  are  shown  in  Figs.  8  and  9  where  y^/kgC  is  plotted  versus 

A. 

k/kg  for  fixed  values  of  o^/ckg  and  y,  and  different  choices  of  the  dimensionless 

2  2  2 

parameters  Apz/pQ  and  ojp/c  kQ.  Note  from  Figs.  8  and  9  that  the  instability 

bandwidth  AK/k^  increases  by  almost  an  order-of-magnitude  as  Apz/pQ  is  increased 
-3  -2 

from  2.8x10  to  2 x lo  .  Also,  the  maximum  growth  rate  increases  substantially 
as  C0p/c2kg  is  increased  from  3.62x10“^  to  1.28  xio-2.  [See  also  Eq.  (2).] 


-13- 


III.  CONTRACT  SUPPORTED  FEL  RESEARCH  PUBLICATIONS/ORAL  PRESENTATIONS 

"Profile  of  a  Relativistic  Electron  Beam  Propagated  through  a  Linear  Wiggler 
and  an  Axial  Guide  Field,"  B.  Hafizi,  G.  L.  Francis,  and  R.  E.  Aamodt,  SAIC 
report  #  SAI C-84- 1428/PR  I -86 ,  to  be  published. 

"Long-Time  Quasilinear  Evolution  of  the  Free  Electron  Laser  Instability  for  a 
Relativistic  Electron  Beam  Propagating  through  a  Helical  Magnetic  Wiggler," 

R.  C.  Davidson  and  Y.  Z.  Yin,  submitted  for  publication  (1984). 

Influence  of  Intense  Equilibrium  Self-Fields  on  the  Spontaneous  Emission  from 
a  Test  Electron  in  a  Realtivistic  Nonneutral  Electron  Beam,"  Ronald  C.  Davidson 
and  Wayne  A.  McMullin,  Phys.  Fluids  27,  1268  (1984). 

"Influence  of  Finite  Radial  Geometry  on  Coherent  Radiation  Generation  by  a 
Relativistic  Electron  Beam  in  a  Longitudinal  Magnetic  Wiggler,"  Ronald  C. 
Davidson  and  Yuan-Zhao  Yin,  Phys.  A30,  3078  (1984). 

"Compton  and  Raman  Free  Electron  Laser*  Stability  Properties  for  a  Cold 
Electron  BeamPropagating  through  a  Helical  Magnetic  Wiggler,"  John  A.  Davies, 
Ronald  C.  Davidson,  and  Georqe  L.  Johnston,  J.  Plasma  Phys.  32_,  in  press  (1985). 

"Cross-Field  Free  Electron  Laser  Instability  for  a  Tenuous  Electron  Beam," 
Ronald  C.  Davidson,  Wayne  A.  McMullin,  and  Kang  Tsang,  Phys.  Fluids  27, 


"Nonlinear  Orbits  in  Free-Electron  Lasers  with  a  Linear  Magnetic  Wiggler  and 
a  Strong  Axial  Magnetic  Guide  Field,"  B.  Hafizi  and  R.  E.  Aamodt,  Phys.  Rev.  A29 
2656  (1984).  - 

"Nonlinear  Traveling-Wave  Equilibria  for  Free-Electron-Laser  Applications," 
Barton  Lane  and  Ronald  C.  Davidson,  Phys.  Rev.  A27,  2008  (1983). 


"Effects  of  Self-Electric  and  Self-Magnetic  Fields  on  the  Raman  FEL  Instability  " 
R.  C.  Davidson,  G.  L.  Johnston,  and  W.  A.  McMullin,  Bull.  APS  28,  1061  (1983). 

Studies  of  a  Tenuous  Beam  Cross-Field  Free-Electron,  R.  D.  Estes,  W.  A  McMullin 
R.  C.  Davidson,  K.  Tsang,  and  A.  T.  Drobot,  Bull.  APS  28,  1142  (1983). 

"Nonlinear  Electron  Orbits  in  Free  Electron  Lasers  with  Axial  Guide  Field  " 

B.  Hafizi  and  R.  E.  Aamodt,  Bull.  APS  28^  1062  (1983). 


"Proposal  for  Continued  Support  of  Research  on  Conversion  Efficiencies  of 
Charged  Particle  Beam  Driven  Free  Electron  Lasers,"  SAI  Proposal  1-254-71- 
840-03,  10  June  1983. 

"Intense  Free  Electron  Laser  Harmonic  Generation  in  a  Longitudinal  Magnetic 
Wiggler,"  Ronald  C.  Davidson  and  Wayne  A.  McMullin,  Phys.  Fluids  26,  840  (1983). 

"Free  Electron  Laser  Instability  for  a  Relativistic  Solid  Electron  Beam  in  a 
Helical  Wiggler  Field,"  Han  S.  Uhm  and  Ronald  C.  Davidson,  Phys.  Fluids  26, 

288  (1983).  — 


-14- 


"Electron  Orbits  in  a  Free-Electron  Laser  with  a  Longitudinal  Magnetic  Wiggler," 
R.  E.  Aamodt,  Phys.  Rev.  A28,  2895  (1983). 

"Nonlinear  Bernstein-Greene-Kruskal  Wave  Equilibria  Subject  to  Global  Energy 
and  Momentum  Conservation  Constraints,"  R.  C.  Davidson  and  K.  T.  Tsang,  J.  Appl . 
Phys.  54,  6284  (1983). 

"Low  Gain  Free  Electron  Laser  near  Cyclotron  Resonance,"  Wayne  A.  McMullin  and 
Ronald  C.  Davidson,  Phys.  Rev.  A25,  3130  (1982). 


"Stability  Properties  of  an  Intense  Relativistic  Non-Neutral  Electron  Ring  in 
a  Modified  Betatron  Accelerator,"  Ronald  C.  Davidson  and  Han  S.  Uhm,  Phys.  Fluids 
25,  2089  (1982). 

"Relativistic  Electron  Beam  Equilibria  with  Non-Circular  Cross  Section,"  Ronald 
C.  Davidson,  Particle  Accelerators  12,  55  (1982). 


"Helically  Distorted  Relativistic  Electron  Beam  Equilibria  for  Free  Electron 
?QinrnQo?\Cati°nS’"  f?ona"ld  c-  Davidson  and  Han  S.  Uhm,  J.  Appl.  Phys.  53 


"Higher  Harmonic  Emission  by  a  Relativistic  Electron  Beam  in  a  Longitudinal 
Magnetic  Wiggler,"  Ronald  C.  Davidson  and  Wayne  A.  McMullin,  Phys.  Rev.  A26 
1997  (1982).  - 

"Proposal  for  Continued  Support  of  Research  on  Conversion  Efficiencies  of 
Charged  Particle  Beam  Driven  Free  Electron  Lasers,  SAI  Proposal  1-254-71- 
830-03,  11  May  1982. 


Stochastic  Particle  Instability  for  Electron  Moti 
Wiggler,  Radiation,  and  Longitudinal  Wave  Fields  " 
Wayne  A.  McMullin,  Phys.  Rev.  A26,  410  (1982) 


on  in  Combined  Helical 
Ronald  C.  Davidson  and 


Free  Electron  Laser  Instability  for 


234M^",er  S  * 


i  n 
24, 


"Theory  of  Free  Electron  Laser  Instability  in  a  Relat 
Beam,  Han  S.  Uhm  and  R.  C.  Davidson,  Phys.  Fluids  24 


vistic  Annular  Electron 
1541  (1981). 


"Proposal  for  Continued  Support  of  Research  on  Conversion  Efficiencies  of 
Charged  Particle  Beam  Driven  Free  Electron  Lasers,"  SAI  Proposal  1-254-71- 
820-01,  20  February  1981. 


Self-Consistent  Vlasov  Description  of  the  Free  Electron  Laser  Instability 
Ronald  C.  Davidson  and  Han  S.  Uhm,  Phys.  Fluids  _23,  2076  (1980). 


-15- 


"Low-Frequency  Flute  Instabilities  in  Intense  Nonneutral  Electron  and  Ion 
Beams,"  Han  S.  Uhm  and  Ronald  C.  Davidson,  Phys.  Fluids  23,  1586  (1980). 

"Proposal  for  Support  of  Research  on  Conversion  Efficiencies  of  Charged 
Particle  Beam  Driven  Free  Electron  Lasers,"  SAI  Proposal  1-254-71-810-03, 

14  May  1980. 

"Semiannual  Status  Report  of  ONR  Contract  No.  N00014-79-C-0555  with  Science 
Applications,  Inc.,  titled:  Conversion  Efficiencies  of  Charged  Particle  Beam 
Driven  Free  Electron  Lasers,"  December  1979. 

"Influence  of  Finite  Radial  Geometry  on  the  Free  Electron  Laser  Instability 
R.  C.  Davidson,  H.  S.  Uhm,  and  R.  Aamodt,  Bull.  APS  24,  1066  (1979). 


Oral  Presentation  to  NRL,  Washington,  by  Richard  E.  Aamodt,  May  1982. 


-16- 


REFERENCES 

1.  R.  C.  Davidson  and  H.  S.  Uhm,  J.  Appl .  Phys.  53,  2910  (1982). 

2.  B.  Hafizi  and  R.  E.  Aamodt,  Phys.  Rev.  A29,  2656  (1984). 

3.  "Influence  of  Finite  Radial  Geometry  on  Coherent  Radiation  Generation  by 
a  Relativistic  Electron  Beam  in  a  Longitudinal  Magnetic  Wiggler," 

R.  C.  Davidson  and  Y.  Z.  Yin,  submitted  for  publication  (1984). 

4.  "Compton  and  Raman  Free  Electron  Laser  Stability  Properties  for  a  Cold 
Electron  Beam  Propagating  through  a  Helical  Magnetic  Wiggler,"  J.  A.  Davies, 
R.  C.  Davidson  and  G.  L.  Johnston,  submitted  for  publication  (1984). 

5.  R.  C.  Davidson  and  H.  S.  Uhm,  Phys.  Fluids  23,  2076  (1980). 


-17- 


FIGURE  CAPTIONS 

Fig.  1.  Normalized  density  D  as  a  function  of  the  distance  z  (in  cm) 

measured  from  the  anode  along  the  drift  tube.  Wiggler  wavelength 
2ir/k  =  3  cm,  y  =  1.68,  guide  field  Hq  =  2.2  kOe,  z  =  <5H/H  =  1/1024, 

resonance  harmonic  number  l  =  2,  radius  of  anode  RQ  =  1.2  cm. 

_2 

Observation  point  for  lower  curve  has  coordinates  (x,y)  =  (0,10  )  cm, 

for  upper  curve  (x,y)  =  (0,0.5)  cm. 

Fig.  2.  Relationship  between  anode  and  typical  gyro-orbits  at  observation 
point  0  with  coordinates  (x,y)  =  (0,0.5)  cm. 

Fig.  3.  Same  parameters  as  Fig.  1,  but  with  e  =  1/128. 

Fig.  4.  Same  parameters  as  Fig.  1,  but  with  e  =  1/64. 

2  2 

Fig.  5.  Normalized  parallel  temperature  T  =  ( ck/2, | } )  (yTz/mc  )  as  a 
function  of  distance  along  the  drift  tube.  Here  z  -  1/128  and 
the  observation  point  has  coordinates  (0,0.5)  cm. 

Fig.  6.  Same  parameters  as  Fig.  1,  but  with  y  =  2.26,  l  =  3,  e  =  1/128.  The 
observation  Doint  for  this  curve  has  coordinates  (x,y)  =  (0,0.5)  cm. 

Fig.  7.  Geometry  and  coordinate  system  for  longitudinal  wiggler  free 
electron  laser. 

Fig.  8.  Plot  of  normalized  growth  rate  Yk/I<QC  versus  k/k^  [Eq.  (2)] 

for  helical  wiggler  FEL  with  eB/mc^k^  =  0.718,  =  3.62x10"^ 

y  =  47.1,  Pq  =  1.286  *10 g-sem/sec,  and  AP2/Pq  =  2.8*  10"^ 

Fig.  9.  Plot  of  normalized  growth  rate  y^/kgC  versus  k/kg  [Eq.  (2)] 

for  helical  wiggler  FEL  with  eB/mc^kg  =  0.718,  co^/c^kg  =  1.28*10”^, 

P  ^ 

Y  =  47.1,  Pq  =  1.286  x  10"^  g-cm/sec,  and  Apz/pg  =  2xlo"^. 
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Fig.  8 
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Fig.  9 


1)2179%  ' 


